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We consider the azimuthal and Pt dependence of hadrons produced in unpolarized Semi-Inclusive 
Deep Inelastic Scattering (SIDIS) processes, within the factorized QCD parton model. It is shown 
\^ , that at small Pt values, Pt < 1 GeV/c, lowest order contributions, coupled to unintegrated (Trans- 

' verse Momentum Dependent) quark distribution and fragmentation functions, describe all data. At 

' larger Pt values, Pt > 1 GeV/c, the usual pQCD higher order collinear contributions dominate. 

Having explained the full Pt range of available data, we give new detailed predictions concern- 
pH ■ ing the azimuthal and Pt dependence of hadrons which could be measured in ongoing or planned 

^ ' experiments by HERMES, COMPASS and JLab collaborations. 



OO 

(N 



PACS numbers: 13.88.+e, 13.60.-r, 13.15.+g, 13.85.Ni 



I. INTRODUCTION 

> . 

' In Ref. Q a comprehensive analysis of Semi-Inclusive Deep Inelastic Scattering (SIDIS) processes within a factorized 
, QCD parton model at 0{a'^) was performed in a kinematical scheme in which the intrinsic tranverse momenta of the 
\^ ' quarks inside the initial proton {k±) and of the final detected hadron with respect to the fragmenting quark (p^) 
I were fully taken into account. The dependence of the unpolarized cross section on the azimuthal angle (j)h between 
' the leptonic and the hadron production plane (Calm effect 0) was compared to the available experimental data, 
' and used to estimate the average values of (fc^) and (p^). These values were adopted in modeling the intrinsic 
(-H , motion dependence of the quark distribution and fragmentation functions. This allowed a consistent description of 
O i' the azimuthal dependence observed by HERMES and COMPASS collaborations in SIDIS off transversely polarized 
' , protons 0, B , with the subsequent extraction P, 0| of the Sivers distribution function ^ . 

In Ref. [l| the main emphasis, following the original idea of Cahn, was on the role of the parton intrinsic motion, 
(-H [ with the use of unintegrated quark distribution and fragmentation functions. That applies to large processes, in a 
• • • kinematical regime in which Pt — Aqcd — k±, where Pt = \Pt\ is the magnitude of the final hadron transverse mo- 
. ^ mcntum. In this region QCD factorization with unintegrated distributions holds Q and lowest order QED elementary 
. processes, £q — *■ £q, are dominating: the soft Pt of the detected hadron is mainly originating from quark intrinsic 
$^ ' motion Ia.l^ ll0l|. ra ther than from higher order pQCD interactions, which, instead, would dominantly produce large 
5^ ; Pt hadrons [n[|llllllli- 

Indeed, a look at the results of Ref. (see, in particular, Figs. 5 and 6) immediately shows that, while the 
inclusion of intrinsic k± and p± - coupled to lowest order partonic interactions - leads to an excellent agreement with 
the data for small values of the transverse momentum Pt of the final hadron, it badly fails at higher Pt'- the turning 
point is around Pt ~ 1 GeV/c. A similar conclusion was drawn in Ref. The large Pt region has been discussed 
at length in the literature and is related to contributions from higher order QCD processes, like hard gluonic radiation 
and elementary scatterings initiated by gluons: these cannot be neglected when Pt S> Aqcd [ill llSll . 

In this paper we start by showing that a complete agreement with data in the full range of Pt can be achieved; for 
Pt < 1 GeV/c we follow the approach of Ref. Q - Pt originated by the intrinsic k± and p± with 0{a^) partonic 
interaction - while in the range of Pt > 1 GcV/c we add the pQCD contributions - collinear partonic configurations 
with higher order [up to ©(a^)] partonic interactions which generate the large Pt- We shall see that indeed most 
available data can be explained; the intrinsic k± contributions work well at small Pt < 1 GeV/c and fail above that, 
while the higher order pQCD collinear contributions explain well the large Pt > 1 GcV/c data and fail, or arc not 
even applicable, below that. The two contributions match in the overlapping region, Pt ^ 1 GeV/c, where it might 
be difficult to disentangle one from the other, as they describe the same physics. Infact parton intrinsic motions 
originate not only from confinement, but also from soft gluon emission, which, due to QCD helicity conservation, 
cannot be strictly collinear. Similar studies, concerning single transverse-spin asymmetries in Drell-Yan and SIDIS 
processes, with separate contributions - TMD quark distributions and higher-twist quark gluon correlations - in 
separate regions, have recently been pubhshed [ig.[l7|. 

Having achieved such a complete understanding of the Pt dependence of the SIDIS cross sections we obtain a full 
confidence on the regions of applicability of the two approaches. We re-analyse the azimuthal dependence of the 
unpolarized cross section - the Cahn effect, described in Ref. - which depends on quantities integrated over Pt- 
The actual data are dominated by the low Pt contributions, and the results previously obtained remain valid; we 
obtain slightly different values of the parameters {k'j_) and (p^). We then consider running experiments (HERMES, 
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COMPASS and experiments at JLab) and physical observables which are being or wiU soon be measured. They are 
mainly in the small Pt regions and we give full sets of predictions for them. 

The plan of the paper is the following: in Section ^ we give a short overview of the kinematics and a collection 
of the basic formulae needed for the computation of the SIDIS cross sections, both in the low Px approach of Ref. 
and in the pQCD large Pt region; in Section Hill we discuss and compare our results for da/dPT, da/d(j)h and 
{cos (ph) with the existing experimental data, over a very wide range of Pt values; in Section llVI we give predictions 
for the forthcoming measurements of da/dPx and {cos4>h) at HERMES, Compass and JLab. Some considerations on 
{cos2(j)h) are made. Finally, in Section we draw our conclusions. 



II. KINEMATICS AND CROSS SECTIONS 



We consider SIDIS processes £p ^ £hX in the 7*p cm. frame, as shown in Fig. ^ The photon and the proton 
collide along the z axis with momenta q and P respectively; the leptonic plane coincides with the x-z plane. We 
adopt the usual SIDIS variables (neglecting all masses): 



s^{P + £y {P + qf = W^ = 2lq2 g2 
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The SIDIS differential cross section can schematically be written in terms of a perturbative expansion in orders of 
as as follows 

da = a° do-Q + a], dai + d(J2 + ■■■ , (2) 

, rfS^^P^ttx ^5^ep-.ehx 

where da is a short hand notation to indicate — ; — ; ,^ „ = x„.s 



dXg. dy dzh c?^Pt dXg. dQ-^ dzh d'^Px 

The first term in this expansion is the lowest order one; in the elementary interaction, £ q —i- £ q, a virtual photon 
with four- momentum q strikes a quark which carries a transverse momentum fcj^ = fcj^(cos if, sin (ys, 0) in addition to a 
fraction x of the light-cone proton momentum. The final detected hadron h originates from the fragmentation of the 
outgoing quark: pj_ is the transverse momentum of h with respect to the direction of the fragmenting quark and z is 
the fraction of the light-cone quark momentum carried by the resulting hadron. Consequently, the detected hadron 
can have a transverse momentum Pt with a magnitude Pt — {k±_) ~ {p±)- Indeed, this is the main source of hadrons 
with a small value of Pt [1 H Elll ■ 
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Such a mechanism translates into a factorized 01 expression for the SIDIS cross section, vahd at all orders in 
(fcx/Q): 



as explained in Ref. 0, where the exact relationships between x,z,pj_ and the observables x^.^z^^Pt are given. 
Notice that at 0{k±/Q) one has x^^- = x, Zh = z and Pt = zk± + Pj_- fq{x,k±) and are the parton 

density and the fragmentation fmiction respectively, for which we assume the usual a;, k± or z,p±^ factorization, with 
a gaussian k± and p± dependence: 

fc^) = f,{x) ^ e-'=i/<'=i> D'^iz,p^) = D'^iz) e-^'i/^^'i) , (4) 

so that 

I d'k^ .fq{x,k^) = f,{x) I d^p^ D';{z,p^) = D'^iz) . (5) 

The integration over d^k± in Eq. induces a dependence on cos<f>h [at 0{Pt/Q)] and on cos2(/)/i [at 0{Pt /Q)"^], 
where <j>h is the azimuthal angle of Pt- The explicit expression, at 0{Pt/Q), is given by: 



dx^^dQ^dzh d^Pr " ^ fq{x,^)D^{zu) 



1 + (1 - y) 



2 



(2-y)yT^(fci)z^P, 



T 



1 



-Pt/{Pt) 



(6) 



where, (P|) = + z2(fc2^). 

Let us now consider the contributions of order Ug , dai in Eq. (0) . We follow the approach of Ref. |0| . The relevant 
partonic processes, shown in Fig. 12 arc those in which the quark emits a hard gluon or those initiated by gluons: 

i* + q^q + 9 i* + q-^9 + q i*+9^q + q- (7) 

It is clear that now, contrary to the lowest order QED process, j* + q q, the final parton can have a large transverse 
momentum, even starting from a coUinear configuration. Such a contribution certainly dominates the production of 
hadrons with large Pt values. 

One introduces the parton variables x' and z' , defined similarly to the hadronic variables x^^ and Zh, 

I -^Bi I k ■ k' Zh 

where k and k' are the four-momenta of the incident and fragmenting partons respectively. ^ and C are the usual 
light-cone momentum fractions, which, in the coUinear configuration with masslcss partons are given by A: = ^P and 
Ph ~ Ck' . We denote by prp (not to be confused with p^) the transverse momentum, with respect to the 7* direction, 
of the final fragmenting parton, Pt = CPt- 

The semi-inclusive DIS cross section, in the QCD parton model with coUinear configuration, can be written, in 
general, as: 



dXg.dydzhd'^P^^ 



- = E / '^^Pt dCS {x,^ - ^x') 6 {zh - Qz') 5^ (P^ - CPt) 



da. 

dx' dy dz' d^PT 

To first order in the partonic cross section is given by [T3 . 
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dcTj 



dx'dydz'(PpT 167r2Q4 



(10) 



where ij denote the initial and fragmenting partons, ij = qq, qg, gq. Inserting the above expression into Eq. yields, 
for the 0{as) cross section: 



2 „2 



dx^Aydzh d^Pr IGn^ x'P} + zl{l - x')Q^ 



dx' 



z/,(l-x')Q2 



with HI 



(11) 



647r 2 {I ■ + jl' ■ k'f + (r ■ kf + (/ ■ k'f 
2, ^ (fc • fc")(fc' • k") 



(i-x')(i-z') + 



1 + jx'z')^ 
(1 -a;')(l - 2') 



+ Sx'z' (1 -y) 



+ 4x'z'(l -?/) cos 20,, I , 



(12) 



T, Mf^" 



6471^2 • kf + [V ■ k" f + {V ■ kf + {I ■ k"f 



— Q 

3 ^ 



(fc • k'){k'k") 



^ (l-x')z' 



+ 8a;'(l-y)(l-/) 



+ 4 ^^^ilJi-^ (2 - y) [x'(l - z') + (1 - x')^'] cos<^. 
+ 4x'(l-2/)(l-z') cos 20;, } , 



(13) 



647r 2 (/ ■ k"f + (/' ■ k'f + (r ■ k"f + (^ ■ k'f 
~Y'^ (fc • fc')(A: • k") 



87ra„ 



[1 + (1 - y)2] [x'^ + (1 - x') 



,^2l 2'V(l-z' 



l\2 



z'{l - z') 



+ 16a;'(l-x-')(l-y) 



^^^^V__^^%^__JJ)_ (2_y) (l_2x')(l-2z') COS0;, 



z'(l - z') 
+ 8x'(l-x')(l-2/) cos 20,, } , 



(14) 
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where we have exphcitely written the scalar products in terms of x\ y, z' and cos 4>h- Notice the appearance of the 
cos (j)h and cos 2(l)h terms: (ph is the azimuthal angle of the fragmenting partons, which, in a collinear configuration, 
concidcs with the azimuthal angle of the detected final hadron. The above expressions agree with results previously 
obtained in the literature [T^ . 

Large values of Pt cannot be generated by the modest amount of intrinsic motion Q ; we expect that Eq. will 
dominantly describe the cross sections for the Icpto-production of hadrons with Pt values above 1 (GcV/c). 

Let us finally briefly consider the contributions of order a1 (NLO), da2 in Eq. ||2J). These, for the production of 
large Pt hadrons, have been recently computed 0,0, resulting in large corrections to the 0{as) (LO) results, and 
leading to a good agreement with experimental data. It would be unnecessarily complicated, for our purposes, to 
take exactly into account these contributions, as we have done for the LO ones, via Eqs. (|ll|) - (|14|l . The most simple 
way of inserting the NLO results in our study is via the K factor, defined as the NLO to LO ratio of the SIDIS cross 
sections. A close examination of the K factor shows a clear dependence on Pt and the other kinematical variables: 
for example, it may be larger than 10 at low Pt and Q^, while approaching unity at larger values (see, for example, 
Fig. 3 of Ref. However, we shall need to use the K factor only in limited ranges of Pt (up to 3 GeV/c at most) 

and (5^, depending on the sets of experimental data we consider: in these limited kinematical regions we find that a 
satisfactory description of the SIDIS experimental data can be achieved by using constant values of A', which will be 
indicated in each case. 

Therefore, we will effectively include 0{a1) contributions in our computations by writing Eq. as 

da = aldcTQ + K dui + ... , (15) 

where da^ and dai are calculated according to Eqs. Q and (|11() respectively. The first term dominates at Pt <1 
GeV/c, and the second one at Pt > 1 GeV/c. 



III. NUMERICAL RESULTS 

In Ref. several sets of experimental data, showing the explicit dependence of the SIDIS unpolarized cross sections 
on the azimuthal angle (f)h and on Pt, were considered. A comprehensive fit, based on Eq. (O ~ or its simplified 
version valid up to 0{k±/Q), Eq. © - was performed, in order to determine the values of (fc^) and (p^), obtaining 

(kl) = 0.25 (GcV)2 {pD ^ 0.20 (GcV)^ . (16) 

These values were assumed to be constant and flavour independent. 

The study of Ref. also showed clearly that Eq. © - zeroth order pQCD with TMD distribution and fragmen- 
tation function - works very well in the Pt — Aqcd — k± region, but fails at larger Pt values, where higher order 
pQCD contributions, with collinear partonic configurations, are expected to take over and explain the data. The 
transition point is around Pt — 1 GeV/c. 

We have redone the analysis of Ref. [2| taking into account, in the appropriate kinematical regions, also the pQCD 
contributions. It turns out that, while a complete description of the data in the full Pt range is possible, a little 
variation is required to the values given in Eq. (I16|l . Actually, the resulting change is included within a 20% variation 
of the parameters, already considered in Ref. 

Figs. l3l7l show our results obtained by adding, according to Eq. H15|) , the contributions of Eq. Q to the contributions 
(computed above Pt = 1 GeV/c) of Eq. ((TT1) . We have used 

(fci) = 0.28 (GcV)2 {pD = 0.25 (GeV)^ , (17) 

again constant and flavour independent. The K factor was fixed to be a constant, with different values according to 
the different Pt and ranges involved. Here and throughout the paper we have adopted the MRSTOl NLO set 
of distribution functions and the fragmentation functions by Kretzer [l9j at NLO. 

Let us comment in greater detail on each single plot. In Fig. Owe compare our results to the EMC measurements 
of the SIDIS P^ distributions (normalized to the integrated DIS cross section) defined as 

I da 1 /" ^^^ip^ihx 

-Jm = o / '^^hdx dydzh- — — — ^5^, (18) 

o'dis ^^ms J dXg. dydzhd^PT 

where the integration covers the x^., y, zt and Pt regions consistent with the experimental cuts: 



Q2 > 5 {GcY/cf Eh>5 (GeV) 0.1 < z,, < 1 0.2 <y< 0.8 



(19) 
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FIG. 3: The normalized cross section da/dP^: the dashed line reproduces the 0{cil) contributions, computed by taking into 
account the partonic transverse intrinsic motion at all orders in {k±/Q), Eq. 0. The solid line corresponds to coUinear and 
pQCD contributions, computed at LO, with a K factor {K = 6) to account for NLO effects, Eqs. Illl l- lll5| l. The data are from 
EMC collaboration measurements {k±) and {p±) axe fixed as in Eq. 1171 . 



The dashed lines reproduce the 0{a^) contribution, computed by taking into account the partonic transverse intrinsic 
motion at all orders in {k±/Q), whereas the solid lines correspond to the SIDIS cross section as obtained by including 
LO corrections and the K factor {K = 6) to account for NLO effects. The two contributions together give a very 
good complementary description of the data over the full Pt domain. Unavoidably, there is a slight mismatch 
at the transition point, Pt = 1 GeV/c, where both contributions somewhat describe the same physics, and some 
kind of average should be performed to avoid double counting. The value if = 6 is the simplest, although rough, 
approximation, in the kinematical range of the data considered here, to the computed K factor (see, Fig. 3 of Ref. 
CTj^js is evaluated starting from Eq. (17) of Ref. [H; pQCD corrections for this integrated quantity are neghgible. 
A similar very good agreement, shown in Fig. 01 is obtained when comparing our computations with the experimental 
measurements of the ZEUS collaboration at DESY Here the SIDIS differential cross section, 

cr^is dPr cTdis J dxg. dQ^ dzh d^Pr 

is obtained by performing the integrations according to the following experimental conditions 

10 < Q2 < 160 (GeV/c)2 7b <W < 175 (GeV) , (21) 

and the K factor is taken to be 1.5 (as shown in Fig. 8 of Ref. 0|). 

In Fig.|31we compare our results to the EMC measurements of the (not normahzed) (j)h distributions 0, proportional 

to 

d(i)h 

The integration regions arc fixed by the conditions 

XF > 0.1 Pt > 0.2 (GcV/c) 0.2 < y < 0.8 > 4(GcV/c)2 , (23) 



- i'^-'y''-''-'''- dx^ dydz.d^Pr - ^''^ 



where xp — ^P^/W and Pl is the longitudinal momentum of the produced hadron relative to the virtual photon. 
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FIG. 4: The normalized cross section da/dPr- the dashed line reproduces the 0(Qs) contribution, computed by taking into 
account the partonic transverse intrinsic motion at all orders in the {k_i_/Q) expansion, Eq. J^; the solid line corresponds to 
the SIDIS cross section as given by LO contributions and a K factor (K = 1.5) to account for NLO effects, Eqs. 
The data are from ZEUS collaboration measurements j2l| . {k±) and {p±) are fixed as in Eq. 1113. 
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FIG. 5: The cross section da/d4'h'- the solid line is obtained by including all orders in {k±/Q), the LO corrections and i 
factor to account for NLO effects. The data are from EMC measurements 0. {k±) and (px) are fixed as in Eq. 11711 . 



K = 6 



This is a quantity integrated over Pt > 0.2 GeV/c; we have used only dao up to Pp ~ 1 GeV/c and added the 
K dai contributions (with K = 6) above that. We notice, however, that the dominant contributions come from very 
low Pt's, while the pQCD contribution is almost negligible. 

Figs.|n|and[7|show our predictions for the average value of cos (f)h compared to the experimental data from the FNAL 
E665 collaboration (TJ {l^p and fid interactions at 490 GeV) and from the ZEUS collaboration (positron-proton 
collisions at 300 GeV) respectively. Here (cos0/i) is defined as 



(cos (j)h) 



J dx^.dQ^dzhd'^Pr cos0^ d^a 



J dx^^dQ^dzfid'^Pr d^a 
where d^a denotes the fully differential cross section 



d^a = 



dXg. dQ"^ dzh d'^Px 



(24) 



(25) 



For the FNAL E665 data sample the integral over Pt runs from P™* to P™°^ ~ 10 GeV/c and the range of the other 



P^"' (GeV/c) 



FIG. 6: {cos (f>h) as a function of P™': the dashed hne reproduces the 0{a^l) contribution, computed by taking into account the 
partonic transverse intrinsic motion at aU orders in {k±/Q); the sohd hne corresponds to the SIDIS cross section as obtained 
by including LO corrections and a K = 6 factor to account for NLO effects. The data are from the E665 coUaboration 



variables is fixed by the following experimental cuts: 



whereas for the ZEUS data sample the integral over Pt runs from P™* to P™°^ ~ 10 GeV/c, within the ranges 



As in the previous case, we have added perturbative corrections only from Pt = 1 (GeV/c), leaving daQ to be the 
only contributing term for values of Pt below 1 (GcV/c). 

The results we obtain are in good qualitative agreement with the FNAL E665 experimental data. As expected, they 
show that the pQCD contributions are very small at low P™* values, but quickly increase as P^"* raises, significantly 
correcting the fast fall of the 0{a'^) term, as shown in Fig. 

Instead, our results disagree with the ZEUS data, especially in the lower range of P^"'. This is surprising and 
would deserve further experimental studies. The cos 4>h modulation, at small Pt values, is a kinematical higher-twist 
effect, and decreases like Pt/Q for growing values of Q, as shown in Eq. Therefore we expect, and indeed we 
find, {cos (f>h) to be much smaller for ZEUS data, which correspond to huge values of Q^, than for E665 results, which 
correspond to much lower values. 

We have also computed, with the same procedure, (cos(20;,)). We have seen that such a dependence can arise 
both, at 0{Pt/Q)^, from intrinsic motion, and, at 0{as), from pQCD corrections. However, there is another non 
perturbative, leading- twist, small Pt source of the cos{2(j)h) dependence, related to the combined action of the Boer- 
Mulders and Collins [l^l effects; this has been recently studied in Ref [2^, where both the kinematical 
contribution and the Boer-Mulders 18) Collins one were studied and found to be of comparable size. Therefore, our 
results, which ignore the Boer-Mulders (X) Collins contribution, can be considered reliable only in the large Pt region; 
indeed, in this region, we find agreement with the ZEUS experimental data of Ref. |22j . as shown in Fig. [S] 



Q2 > 3 (GcV)2 , 300 <W^ < 900 (GcV)^ , 

60 < J/ < 500 (GeV) , Eh > 8 (GeV) , O.K y < 0.85 



(26) 



0.01 < x^^ < 0.1 , 0.2 <y < 0.8 , 

0.2 <zh< 1.0 , 180 < 0^ < 7220 (GcV)^ 



(27) 
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FIG. 7: {cos (fih) as a function of P™': the dashed hne reproduces the 0{a^l) contribution, computed by taking into account the 
partonic transverse intrinsic motion at aU orders in {k±/Q); the sohd hne corresponds to the SIDIS cross section as obtained 
by including LO corrections and a K — 1.5 factor to account for NLO effects. The data are from the ZEUS collaboration |2^. 




P^"' (GeV/c) 



FIG. 8: {cos{2(j>h)) as a function of P™* as obtained by including LO corrections and a. K = 1.5 factor to account for NLO 
effects. The data are from the ZEUS collaboration |23 . 

IV. PREDICTIONS FOR FORTHCOMING MEASUREMENTS 

New data are expected from ongoing measurements or data analysis at HERMES, COMPASS and JLab. They 
coneern dominantly the small Pt (and, hopefully, large enough Q^) region, where we have seen that the simple partonic 
approach, with unintegrated distribution and fragmentation functions, can give a very satisfactory description of the 
available data. We can easily give detailed predictions which can soon be tested, allowing a further check on the role 
of intrinsic motions in affecting physical observables. We consider the SIDIS cross sections and the average value of 
cos(f>h- Also {cos2(j}h), keeping in mind the comments at the end of the previous Section, is computed. 

In Fig. El we plot the SIDIS cross section, for tt"*" production at HERMES, as function of one variable at a time, 
either Zh, x^., y or Pt; the integration over the unobserved variables has been performed consistently with the setup 
of the HERMES experiment, which studies the scattering of positrons at piah ~ 27.57 GeV/c against a fixed hydrogen 
gas target: 

Q2 > 1 {GeV/cf > 10 GeV^ Pt > 0.05 GeV/c 

0.023 < x^^ < 0.4 0.2 < zh < 0.7 0.1 < y < 0.85 2 < Eh < 15 GeV . (28) 

In these kinematical regions the cross section is heavily dominated by the dao term of Eq. (O , computed according 
to Eq. (PJ; (JDis is computed according to Eq. (17) of Ref. 0. We also evaluate the average value of cos (j)h in the 
same kinematical region, shown in Fig. IIUI and of cos2(j)h, shown in Fig. 1111 The latter, however, can only be taken 
as a partial (higher-twist) contribution to the real (cos20/i), as explained at the end of the last section. 
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We notice that very similar predictions, for the cross section, {cos (ph) and (cos2(/)/j), are obtained for 7r° and tt 
production, which we do not show. 

The COMPASS experiment at CERN collects data in ^id ^h^X processes at piab — 160 GeV/c, covering the 
following kinematical regions: 

> 1 (GeV/c)2 > 25 GeV^ Pt > OA GeV/c 
E,, < 15 GeV 0.2 < < 1 O.K y < 0.9 . (29) 

Fig. 1121 shows our corresponding predictions for the SIDIS cross section - for the production of positively charged 
hadrons - as a function of the kinematical variables Zh, x^., y and Pt, as obtained from Eq. Q. Notice that we 
neglect nuclear corrections and use the isospin symmetry in order to obtain the parton distribution functions of the 
deuterium. Similarly to the HERMES case, the perturbative QCD corrections are negligible. Predictions for {cos (j)h) 
are presented in Fig. Upland the higher-twist contributions to (cos20?i) in Fig.^] Similar results hold for negatively 
charged hadrons. 

Finally, JLab collects and will collect data in the collisions of 6 and 12 GeV electrons from a fixed He'^ target. In 
this case the relevant kinematical regions are the following: 

Q2 > 1 (GeV/c)2 W^>AGcY^ Pt > 0.1 GeV/c 

1.0 <Eh< 3.5 GeV 0.4 < z^ < 0.7 0.1 < x^^ < 0.6 OA < y < 0.85 . (30) 

Our results for the corresponding SIDIS cross sections, {cos (j)h) and (cos20h) are shown in Figs. IT ^ 1161 and 1171 
respectively: as for the previous experiments, also the JLab data are dominated by 0{a^) terms and are almost 
completely insensitive to the LO and NLO perturbative QCD corrections. Again, we show results for tt"*" production, 
but very similar ones hold for tt~ and tt'^, which we do not show. 

All these results depend on intrinsic momenta, both k± in the partonic distributions and in the quark fragmen- 
tation. This is obvious for quantities like da/dPr and {coscfih) which could not even be defined, at 0{a^), without 
intrinsic motion (and pQCD corrections are negligible for the experiments we consider). However, this is also true for 
the differential cross-sections in z, x^. and y: although they get contributions from intrinsic motion only at 0{Pt /QY' , 
as one can explicitly see from Eq. (|^, these contributions can be sizeable in the kinematical domains of HERMES, 
COMPASS and JLab. 



V. CONCLUSIONS 

We have considered the azimuthal and Pt dependence of SIDIS data from low to large Pt; they both cannot 
be explained in the simple parton model (©(a^)) with collinear configurations. They can originate from intrinsic 
motions and/or pQCD corrections. The outcome of our analysis turns out to be very simple: up to Pt < 1 GeV/c 
the simple parton model with unintegrated distribution and fragmentation functions explains the data and leads to 
a good evaluation of (fc^) and while at larger Pt values, Pt > 1 GeV/c, the perturbative QCD contributions 

originating from hard gluonic radiation processes and elementary scattering initiated by gluons, are dominant. 

Having clearly established the complementarity of the two approaches, and in particular having gained full confidence 
in the domain of applicability of the unintegrated parton model, we have given predictions for the cross sections and 
for the values of (cos(/)/j), as they will soon be measured by HERMES, COMPASS and JLab collaborations, mainly 
in the low Pt region. These new data will be a very important tool to test our knowledge of the intrinsic partonic 
internal motion and on the TMD quark distribution and fragmentation functions. 
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FIG. 9: Predictions for the normalized SIDIS cross section corresponding to tiie production of tt""" as it will be measured by 
the HERMES collaboration in the forthcoming future. The solid lines correspond to the SIDIS cross section as obtained by 
including all orders in the {k±/Q) expansion. Notice that QCD corrections have no influence in this range of low Pt's. 
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FIG. 10: Predictions for (cost^h) corresponding to the production of n'^ as it will be measured by the HERMES collaboration 
in the forthcoming future. The solid lines correspond to (cos^h) we find by including all orders in the {k±/Q) expansion. 
Notice that QCD corrections have no influence in this range of low Pt's. 
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FIG. 11: Predictions for (cos 24>h) corresponding to the production of tt"*" as it will be measured by the HERMES collaboration 
in the forthcoming future. The solid lines correspond to (cos2(/)fi) we find by including all orders in the {k±/Q) expansion. 
Notice that QCD corrections have no infiuence in this range of low Pt's. 
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FIG. 12: Predictions for the normalized SIDIS cross section corresponding to the production of positively charged hadrons as 
it will be measured by the COMPASS collaboration in the forthcoming future. The solid lines correspond to the SIDIS cross 
section as obtained by including all orders in the {k±/Q) expansion. Notice that QCD corrections have no influence in this 
range of low Pt's. 
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FIG. 13; Predictions for (cos^h) corresponding to the production of positively charged hadrons as it will be measured by the 
COMPASS collaboration in the forthcoming future. The solid lines correspond to (cos(j!>h) we find by including all orders in 
the {k±/Q) expansion. Notice that QCD corrections have no influence in this range of low Pt's. 
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FIG. 14: Predictions for {cos2(fih} corresponding to the production of positively charged hadrons as it will be measured by the 
COMPASS collaboration in the forthcoming future. The solid lines correspond to (cos20h) we find by including all orders in 
the {k±/Q) expansion. Notice that QCD corrections have no influence in this range of low Pt's. 
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FIG. 15: Predictions for the normalized SIDIS cross section corresponding to the production of vr"*" as it will be measured 
by the JLab collaboration in the forthcoming future. The solid lines correspond to the SIDIS cross section as obtained by 
including all orders in the {kj_/Q) expansion. Notice that QCD corrections have no influence in this range of low Pt's. 
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FIG. 16: Predictions for (cos^h) corresponding to the production of vr^ as it will be measured by the JLab collaboration in 
the forthcoming future. The solid lines correspond to (cos^/i) we find by including all orders in the {kj_/Q) expansion. Notice 
that QCD corrections have no influence in this range of low Pt's. 
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FIG. 17: Predictions for (cos2(^fi) corresponding to the production of 7r+ as it will be measured by the JLab collaboration in 
the forthcoming future. The solid lines correspond to (cos2(/)fi) we find by including all orders in the {k±/Q) expansion. Notice 
that QCD corrections have no influence in this range of low Pt's. 



